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Abstract 



The implementation of QED initial-state radiative corrections in the process 
of four-fermion production at LEP200 and higher-energy e + e~ colliders is 
discussed. Because of the presence of charged-current processes, this is a 
nontrivial problem, and we compare our approach with other existing treat- 
ments. We describe the Monte Carlo algorithm used for the generation of 4- 
fermion events with photon bremsstrahlung. Comparison between our event 
generator and semi-numerical calculations are presented, as well as predic- 
tions for W and Z pair related cross sections. 
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1 Introduction 



In this letter we discuss the implementation of initial-state QED radiation 
(ISR) in the reaction 

e + e~ — > 4 fermions . (1) 

We shall discuss our method, present a variety of results for the energy regime 
of LEP200 and slightly above, and, where possible, make comparisons with 
other existing results. The purpose of this study is to provide an efficient 
and flexible calculational tool, in the form of a Monte Carlo event generator, 
EXCALIBUR, that incorporates the expected dominant physics effects. There- 
fore, this tool should be useful for all kinds of physics studies at LEP200 and 
beyond. We start by outlining our method. 

2 The method 

In a recent paper Q, we studied all electroweak (EW) e + e~ processes leading 
to a final state of four effectively massless fermions, without regard for the 
effects of photon emission. In the case of four quarks, one would like to 
add the concomitant QCD production channel, and also the production of 
a quark pair and two gluons, since both different final states will appear as 
jets. This QCD-extended treatment is now also availabe ||. 

The motivation for our considering 4-fermion final states is that the single 
or pairwise production of vectorbosons, at LEP200 energies or beyond, typ- 
ically manifest themselves as 4-fermion production. A priori, one does not 
know the relative sizes of the production via an intermediate step of vector 
boson production (the 'signal'), or via direct production (the 'background'), 
and, moreover, these two alternatives interfere; for the interpretation of fu- 
ture experiments this knowledge is requisite. The results of Q show that, 
for a number of final states, background effects can be as large as radiative- 
correction effects. 

In a general procedure was developed in which any 4-fermion final 
state can be chosen, and (weighted) Monte Carlo events efficiently generated. 
In this way, every possible distribution or cross section in these processes 
becomes calculable. The only restriction in is that charged particles in the 
final state be, experimentally, visible and distinguishable, i.e. they must be 
produced with sufficient energy, and at nonnegligible angle to the beams or to 
each other. Under this condition, one avoids collinear singularities: therefore, 
the fermions can be treated as massless, which considerably accelerates the 
calculation []. When one considers only collinearly convergent diagrams, we 

1 Note that signals with an intermediate Higgs therefore fall outside of our scope. 
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may even omit our visibility requirements. Amongst others, this case occurs 
for the W pair signal reaction 

e + e~ — ► W + W~ — > 4 fermions . (2) 

In the calculation of the Born diagrams that are needed in the event gen- 
erator, a number of input parameters can be chosen independently, to wit, 
the electroweak couplings parametrized by a (running) a and sin 2 ^, the 
boson masses m w and m Z) and the corresponding total widths Y w and Y Z - 
For instance, one has the freedom to choose cos9 w different from m^/mz, 
and to take arbitrary values for the widths, which may also be assigned an 
energy dependence. This convention, in which a occurs as an overall factor 
in the matrix element, automatically ensures the unitarity cancellations, so 
that the value of a may be changed at will. 

Although the event generator of |]J is suitable for many signal-versus- 
background studies, its usefulness is increased considerably when the most 
important radiative correction effects are incorporated. Let us discuss how 
this can be achieved in a practical way. 



3 Inclusion of radiative corrections 

To start, we focus on the reaction 

e + e~ -> W + W~ , (3) 

with stable TV's (T w = 0), and on the reaction in eq. (H), with decaying Ws 
(IV > 0). The radiative-correction problem is much more involved than for 
the, by now familiar, case (see, e.g. 0) 

e + e- -> Z° -> . (4) 

In the latter reaction, photonic and weak corrections can be separated. Thus, 
the weak one-loop corrections lead to a modification of the Born cross section 
into a 'dressed' Born cross section. To this latter, the sizeable QED initial- 
state corrections can be applied by means of the structure-function method. 
These QED corrections incorporate O (a) and O (a 2 ) leading-log (LL) and 
subleading terms, while the leading higher-order soft-photon contributions 
are implemented by exponentiation. 

Despite many attempts, a simple dressed Born cross section for the re- 
action of eq. (0) has not be found, so that the above procedure for the Z 
cannot be applied: the weak corrections do not decouple from the the QED 
ones in W pair production. A cognate difficulty is that a division of pho- 
ton emission into initial- and final-state radiation is meaningful for reaction 
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Table 1: A summary of references to radiative-correction studies of W pair 
or 4-fermion production. Either electroweak or LL QED ISR have been con- 
sidered. Event-generator approaches are mentioned explicitly. All references 
contain numerical results, except ref. 
most complete treatment. 



14 1 which devises a strategy for the 



(H) but not for (y), since the two sets of Feynman diagrams are not sepa- 
rately gauge-invariant. A priori, then, it is meaningless to consider those 
for initial-state radiation in isolation. An elegant trick has been proposed in 
0] to introduce a gauge-invariant definition of initial- and final-state terms, 
based on adding and subtracting extra terms in the radiative matrix ele- 
ment. Although solving, in a sense, the problem of gauge invariance, there 
is as yet no proof that the initial-state radiation terms thus obtained yield a 
quantitatively good description. Since the LL terms in the cross section are 
anyhow gauge-invariant, the method of [|J is correct for these terms, but the 
subleading ISR terms now contain some arbitrariness. 

In order to facilitate the discussion that follows, we find it useful to sum- 
marize in table 1 what has been done in the literature, and what some authors 
hope to achieve. We also indicate the position of the present paper in this 
welter of possibilities. 

For stable W's, complete one-loop EW effects have been calculated in a 
number of papers || || [7|, ||. The results of refs. [^] and || are in perfect 
agreement. To obtain the complete O (a) correction, the effect of emission 
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of a single hard photon has to be included |§. An event generator based on 
this calculation also exists [jn| , whereas recently another event generator has 
been constructed where, for collinear photons, O (a 2 ) QED corrections are 



also considered 11 



In order to asses the importance of higher-order photonic effects, a LL 
calculation up to 0(a 2 ), including exponentiated soft-photon effects, was 
already carried out in an early for the total cross section and for an angu- 
lar distribution ||12|| . Somewhat later, the O (a) term of such a QED LL 



treatment was replaced by a full O (a) EW calculation [^, IB]. The differ- 
ence between these two approaches amounts to less than 1% in the LEP200 
energy range. 

For decaying W's, the radiative corrections are far less known. The full 
O (a) EW corrections have not yet been computed. Questions of strategy, 
and theoretical issues like gauge invariance, have been addressed in the liter- 
ature, which may form the basis for a conclusive calculation ]14[|. What has 



been done, are ISR calculations to the 4-fermion final state produced by the 
WW signal. In ref. ||, the LL plus exponentiated soft-photon contributions 
are evaluated semi-analytically. For the total cross section, a 3-dimensional 
integral has to be done numerically, after a number of analytic integrals 
have been performed. Thus, only very specific distributions can be obtained. 
The influence of sub leading terms in the ISR are studied, but, since the 
above-mentioned trick has been used in Q in order to give a gauge-invariant 
definition of ISR, the quantitative meaning of precisely these terms is obscure 
until a full calculation, including final-state radiation (and interference) has 
been achieved. Another approach uses LL QED corrections to the signal 



diagrams, leading to an event generator for the three signal diagrams [15 



In the present paper, we want to be able to study any experimental distri- 
bution f] and all background effects, and we want to incorporate the dominant 
radiative correction effects. This leads to an event generator that can han- 
dle all diagrams leading to a specified 4-fermion final state (with, of course, 
the option of a restriction to the signal diagrams), and that incorporates 
the LL O (a) and O (a 2 ) ISR, with exponentiation of the remaining soft- 
photon effects. We shall now discuss how the initial-state radiation effects 
are incorporated in our Monte Carlo event generator 



2 We except the transverse-momentum distribution of the bremsstrahlung. 



4 



4 Generating initial-state radiation 



In order to upgrade the event generator of ref. |]] with QED ISR, the fol- 
lowing description of the radiation is used. Each of the incoming fermions 
is assumed to have its energy degraded by an amount of bremsstrahlung. 
Under the assumption that the bremsstrahlung photons are emitted parallel 
to the radiating beam, the energy distribution of the fermion after radiation 
is described by the 'structure' function 



"8vr 2 



I _|_ 3-£ 2 

log x + 4(1 + x) log(l — x) + 5 + x 



x 



7T 




P = -(L-1) , L = log(^) , (5) 



where x is the fermion's energy in units of the beam energy, m e is the electron 
mass, 7s is Euler's constant, and Q 2 is some appropriate energy scale. This 
is a LL O (a 2 ) structure function with exponentiated soft-photon effects. 
Subleading logarithmic terms are not considered, since one needs for this the 
non leading logarithmic terms from the one loop EW corrections to the Born 
cross section of reaction in eq. (|I|). Another important point is the choice of 
Q 2 . Formally, a change in Q 2 is a subleading effect, but its precise value is 
of course a matter of numerical concern. We use Q 2 = s, the total energy 
squared, which is known to be acceptable [O]. It is a unique advantage of 



the Monte Carlo approach that, if we wish, the Q 2 can even be determined 
on an event-by-event basis. Finally, the structure functions <&(xi) for the 
incoming e + (with original momentum p±, degraded to x±pi) and $(22) f° r 
the incoming e~ (with momentum P2, degraded to X2P2), can be convoluted: 
with 

s' = (xipi + X2P2) 2 = xix 2 s , (6) 

we arrive at the 'flux function' 

1 1 

G{s'/s) = J J dxi dx2 $(#2) 5(x\X2 — s'/s) (7) 



which enables one to write the total radiative cross section as 

1 

a(s) = J dz G(z)a (zs) , (8) 


where <Jq is the nonradiative cross section. Incidentally, the form of G(x) 
is, to given order in a, quite close to that of $(x), where a is replaced by 



5 



2a ]T3|]. In LL approximation the flux function Gd of ref. || is related to 



the choice made in eq. @. We have opted for the use of structure func- 
tions rather than that of a flux function for the following reason. Assuming 
that the four-momentum lost to radiation is lightlike (as one usually does 
in problems such as this one), the total energy loss of the beams must be 
equal, in the flux-function formalism, to (1 — s'/s)a/s/2; for instance, this 
identification is made in M. The energy loss can, however, be quite different 
in the structure-function formalism, where the lost momentum is the sum 
of two lightlike vectors; and this is the more realistic approach, since radia- 
tion from the two beams tends to be contained in two narrow, back-to-back 
cones. Another way to appreciate the difference is to note that, when all the 
lost momentum is lumped into a single lightlike vector, events with small s' 
will always be boosted away from the lab frame, whereas in the structure- 
function formalism they can easily be at rest in the lab frame. 

We are now in a position to decribe the radiation algorithm. To start, 
two values for x\ and X2 are generated, each with a probability distribution 

= ^(l-a;)^- 1 , 
Pi 

7i \ m\ J 

where we have introduced yet another scale Q\. We then compute x\p\ and 
X2P2, so that the total momentum of the 4-fermion system is determined. 
This allows us to move over to the centre-of-mass frame of the 4-fermion 
system; here, the total energy is v? rather than y/s, and, importantly, the 
beam directions are the same as in the lab frame. This is due to the assump- 
tion that the bremsstrahlung is strictly collinear: relaxing this assumption 
would lead to, as yet nearly insuperable, complications in the Monte Carlo. 

We now generate a weighted Monte Carlo 4-fermion event by the pro- 
cedure described in [|T|. The resultant momenta are then boosted back to 
the lab frame. The final event weight is then the product of the '4-fermion 
weight' defined in jjj], and the 'radiation weight' 



$i(xi) $i(x 2 ) 



w rad = ^ /„ ^ /„ n ' ( 10 ) 



As usual, the Monte Carlo cross section and its estimated error are then 
extracted from the mean and the variance of the weight distribution. We 
want to stress that the adoption of a structure function different from that of 
eq. @ is trivial in our approach, since it only entails modifying the definition 
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of w ra( ^: in this way, we can, for instance, perform delicate checks with the 
results of ref. [[§]. 

A few remarks are in order here. In the first place, note that the use 
of two structure functions with terms up to, say 0(a), is not equivalent 
to a calculation based on a flux function to O (a), since the product of the 
two structure functions contains some O (a 2 ) terms. Of course, with the 
above algorithm we may also settle for generating only a single value for 
z = x±X2, and proceeding with the generation of the 4-fermion final state 
using the reduced energy s' = zs. Since, however, the Lorentz boost is then 
not determined, we can only compute the total cross section in that case, and 
no differential ones. Nevertheless, we have applied this procedure in order to 
compare our event generator with the semi-analytical results of ref. 

Secondly, we have to discuss the choice of Q\. In principle, it may be 
chosen at will, since the Monte Carlo cross section does, formally, not depend 
on it. However, it is easily seen that u> rac i will diverge as x\ or x 2 approache 
one, unless Q\<Q 2 . We therefore must, in the generation of events, choose 
Q\ to be the minimum possible value for Q 2 . If one desires to use a fixed Q 2 
scale (as we have done in this paper) this poses no problem, and one simply 
puts Q\ = Q 2 ; but for a study where Q 2 depends on the particular event 
generated, some care must be taken. 

Finally, it should be realized that the various kinematical cuts described 
above must also be boosted to the 4-fermion rest frame. For invariant-mass 
cuts this is no problem, but the energy and angular cuts require some atten- 
tion. Given two values x\ and x 2 , we know the relativisitc velocity of the 
Lorentz boost to the 4-fermion centre-of-mass frame (CMF): 

f3 = ( Xl -x 2 )/(xx +x 2 ) • (11) 

If the scattering angle 9 of a (massless) particle in the lab frame is restricted 
between 

c < cosO < ci , (12) 

we may then compute the bounds in the CMF on the CMF scattering angle 
0': 

Co - P ^ Q , ^ Ci -p 

<cosfl < . 13 

l-co/3 l-ci/3 v ; 

Similarly, suppose its energy E in the lab is restricted by 

E > £ min . (14) 

Its energy E' in the CMF now depends on both the energy and the angle in 
the lab frame, so that a CMF energy cut is complicated; we replace it by its 
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lower bound (the minimum over all scattering angles): 

E' > (l-max(c /3, Cl /3)) . (15) 

Since this cut is somewhat looser, some particles may end up with an energy 
lower than E m ^ n ; this means that an additional number of generated events 
has to be rejected. It must be stressed that these cuts would become un- 
practically complicated if the bremsstrahlung were also to have a transverse 
momentum component. 



5 Results and conclusions 

We shall now discuss a number of results from EXCALIBUR. In the first place, 
we have to establish agreement, where possible, with the results of ref. ||. 
To this end, we must of course make sure to use the same electroweak input 
parameters. In the structure functions and the flux functions we use the 
low-energy value 

«str.f., flux = (137.036)- 1 . (16) 

In the rest of the calculation, we take the input parameters such that impor- 
tant electroweak corrections to the reaction are effectively incorporated. 
In the notation of ref. [13], the amplitude M for reaction (|3]) can be divided 



into two gauge-invariant parts: 

M =2iiik M ' +e2Mo - (17) 

where Mj is the purely V — A part, and Mq a purely vector-like part. The 
number e 2 = Ana we choose such that the running value of a at LEP200 
energies is obtained, and sin 2 6 W we fix by 

a = G^m 2 w 
2sin 2 ^ vrV2 ' 

This gives us the parameters also used in ||. We also adopt the values used 
there for the boson masses and widths. Numerically, we then have 

a = (127.29)" 1 , sin 2 ^ = 0.2325 , = 1.16635 10~ 5 GeV~ 2 , 
m z = 91.173 GeV , T z = 2.4971 GeV , 

m w = 80.220 GeV , T w = 2.033 GeV . (19) 

We use, moreover, boson propagators with energy- dependent widths; their 
denominators have the typical form s — m 2 + isT/m. In [O], a discussion of 
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an 'improved Born approximation' is given. In our approach, where we have 
both signal and background effects, we can retain of this discussion only the 
above-mentioned tuning of a and sin 2 9 W . In particular, the effects of the 



Coulomb singularity in the WW system |L6], |lj], is left out. 

We are now ready to present some saliant results. We have chosen the 
explicit process 

e + e~ — > e~v e ud , (20) 
which contains the WW pair production signal, as well as a nonnegligible 



background. In the first place, we must reproduce the results of j|, [18| where 
we can. These are the Born cross section oq\ the total cross section with 
ISR from ref. [18] and the average 'energy loss', which in Q] is defined as 



^(1 — X1X2), which we denote by e. This is (see above) not the real energy 

loss 2 (2 — x\ —X2), denoted by e; but which of the two quantities is actually 
the most relevant in the measurement of the W mass is, of course, determined 
by the prospective data analysis. 

We have performed a number of different Monte Carlo studies, under 
different strategies, which we now list. 

• WW,f,l: the WW signal diagrams only, with the flux-function ap- 
praoch to 0(a)] no cuts. 

• WW,s,l,a: the WW signal diagrams only, with structure functions 
that are simply the flux function of to O (a), in which 2a is replaced 
by a; no cuts. 

• WW,s,l,b: the WW signal diagrams only, with the structure func- 
tions of eq. (^) where the last (L 2 ) terms are left out; no cuts. 

• WW,f,2: the WW signal diagrams only, with the complete O (a 2 ) flux 
function of ref. || (which, unfortunately, is not given explicitly); no 
cuts. 

• WW,s,2: the WW signal diagrams only, with structure functions as 
given in eq. (|); no cuts. 

• WW, cuts: like the previous case, except that now we also impose the 
following cuts: 

E e -, u ,d > 20 GeV, I cos6 l e - i Uj j| < 0.9, | cos (ud)\ < 0.9, m(ud) > 
10 GeV. 

• all, cuts: like the previous case, except that now also all the background 
Feynman diagrams are taken into account. 
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The various results are given in table 2, where the numbers taken over from 
ref. H, [18| are indicated by the subscript 4. The cross sections are given in 
picobarns, and the energy losses in GeV. 

Similarly in table 3 results are given for Z-pair production. The choice 
of the second energy was determined by the condition that the velocities 
of the Z's would be the same as for the W's at \/s = 190 GeV, when 
there is no radiation present. The cuts used are E( a n particles) > 20 GeV, 
| cos 9 ( a ii P articies)\ < 0.9, m(e+ e -) and m(uu) > 10 GeV and | cos (uu)\ < 0.9. 

In discussing the results, let us first notice that the cross sections from 



the event generator and the semi-analytical approach j|, pl| agree. For ctq 
it is clear when one takes the numbers .60078 and .67930 from fT%| , for a 
it follows from the flux function comparisons in table 1. These comparisons 
agree within 1 standard deviation. For e and 64 there are differences up to 
five standard deviations. The error on the results for q of ref. M is however 



not available |L8|, so that no conclusion can be drawn. 

The structure function method differs slightly from the flux function 
method as comparisons between the first and second rows of table 1 shows. 
The O (a 2 ) results from EXCALIBUR and ref. [§] differ at the 2% level for 
the energy losses, but this may be due to the not specified form of O (a 2 ) 
corrections in Q. 

The inclusion of cuts and the inclusion of more diagrams affect both cross 
sections and energy losses. 

Since the proposed direct reconstruction method for the W mass suffers 
from a shift in M w due to the radiated energy fl9j , a precise knowledge of 
e and e is warranted. Our results show that the precise treatment of ISR, 
the choice of cuts and the neglect of diagrams all affect the energy losses. 
The first effect was also found in ref. M, where also the influence of the 
Coulomb singularity is discussed. The two other effects show that a Monte 
Carlo treatment allowing for cuts and being able to include all diagrams is 
indispensable. 

The results for Z-pair production again show the effects of cuts and the 
inclusion of background diagrams. Although the second energy is tuned in 
a way that comparisons to W-pair production may make sense, the energy 
losses divided by the total energy are different for these two processes. 
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strategy 


(T <T 4 


a 


e 4 


e 


e 


v/i =176 GeV 


WW,f,l 


.60111 
.00032 


.50504 


.50490 
.00032 


1.168 


1.162 
0.002 




WW,s,l,a 


» 




.50484 
.00033 




1.172 
0.002 


1.175 
0.002 


WW,s,l,b 






.50175 
.00098 




1.167 
0.006 


1.170 
0.006 


WW,f,2 


ii 


.50315 




1.200 






WW,s,2 


ii 




.50258 
.00097 




1.178 
0.006 


1.181 
0.006 


WW,cuts 


.44651 
.00092 




.37737 
.00091 




1.192 
0.007 


1.195 
0.007 


all, cuts 


.45011 
.00097 




.37926 
.00095 




1.149 
0.007 


1.152 
0.007 


v/i =190 GeV 


WW,f,l 


.67911 
.00038 


.60733 


.60750 
.00039 


2.108 


2.092 
0.003 




WW,s,l,a 


ii 




.60706 
.00039 




2.111 
0.003 


2.118 
0.003 


WW,s,l,b 


ii 




.60330 
.00117 




2.098 
0.010 


2.106 
0.010 


WW,f,2 


ii 


.60696 




2.171 






WW,s,2 


ii 




.60630 
.00117 




2.124 
0.010 


2.132 
0.010 


WW,cuts 


.48289 
.00107 




.43124 
.00106 




2.188 
0.013 


2.196 
0.013 


all, cuts 


.49316 
.00115 




.44164 
.00114 




2.136 
0.013 


2.144 
0.013 



Table 2: Results on radiatively corrected cross sections and average energy 
losses, under various calculational strategies, for the process e + e~ — > e~v e ud. 
The second line in each entry is the estimated Monte Carlo error. 



13 



strategy 




a 


e 


e 




V 


s = 


190 GeV 






ZZ,s,2 


.70549 10 


-2 


.54632 10" 2 


0.744 


0.745 




.00040 




.00060 


0.002 


0.002 


ZZ,cuts 


.52192 10' 


-2 


.40742 10" 2 


0.732 


0.734 




.00060 




.00065 


0.003 


0.003 


all, cuts 


.92509 10" 


-2 


.78393 10" 2 


1.929 


1.940 




.00182 




.00182 


0.016 


0.016 






=2 


L5.942 GeV 


ZZ,s,2 


.93245 10 


-2 


.83903 10- 2 


2.614 


2.624 




.00044 




.00067 


0.008 


0.008 


ZZ,cuts 


.62890 10- 


-2 


.57262 10" 2 


2.699 


2.709 




.00079 




.00084 


0.011 


0.011 


all, cuts 


.99837 10' 


-2 


.92763 10" 2 


3.503 


3.526 




.00194 




.00198 


0.021 


0.021 



Table 3: Results on radiatively corrected cross sections and average energy 
losses, under various calculational strategies, for the process e + e~ — > e + e~~uu. 
The second line in each entry is the estimated Monte Carlo error. 
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